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Dynamic  Multiscale  Simulation  of  Polyelectrolyte  Nanoassemblies 

ABSTRACT 

The  goal  of  this  project  is  to  develop  a  fundamental  understanding  of  how  fluid  flows,  electric  fields,  and  surface  patterning  affect 
polyelectrolyte  adsorption,  self-assembly,  and  transport  near  solid  surfaces.  Polyelectrolytes  are  charged  polymers  whose  behavior  is 
central  to  Anny-relevant  technologies  (e.g.,  biochemical  sensors,  surface  functionalization)  but  poorly  understood.  Brownian  dynamics 
simulations— a  coarse-grained  simulation  method  suitable  for  describing  non-equilibrium  polymer  behavior— have  been  applied  to  study 
polyelectrolyte  adsorption  in  shear  flow,  polyelectrolyte  adsorption  onto  patterned  surfaces,  dendrimer  and  hyperbranched  polymer 
adsorption,  and  polyelectrolyte  electrophoresis.  Along  with  complementary  continuum-level  models,  the  simulations  yield  insight  into  how 
to  design  fluid  flows,  electric  fields,  and  surface  patterning  to  control  the  behavior  of  single  polyelectrolyte  molecules.  The  knowledge 
gained  lays  the  foundation  for  future  studies  involving  multiple  molecules  and  more  complex  physical  phenomena. 
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Dynamic  Multiscale  Simulation  of  Polyelectrolyte  Nanoassemblies 
Satish  Kumar,  University  of  Minnesota 

Statement  of  Problem  Studied 

Background 

Polyelectrolytes  are  charged  polymers  (e.g.,  DNA  and  proteins)  whose  adsorption,  self- 
assembly,  and  transport  are  poorly  understood  phenomena  relevant  to  a  number  of 
important  technologies  such  as  biochemical  sensors,  nanofluidics  and  microfluidics, 
biocompatible  coatings  and  materials,  separation  membranes,  and  microelectronic  and 
photonic  devices.  The  goal  of  this  project  is  to  develop  a  fundamental  understanding  of 
how  fluid  flows,  electric  fields,  and  surface  patterning  (chemical  and/or  topographical) 
affect  polyelectrolyte  adsorption,  self-assembly,  and  transport  near  solid  surfaces.  Such 
fundamental  understanding  will  enable  the  rational  design  of  fluid  flows,  electric  fields, 
and  surface  patterns  for  creation  of  polyelectrolyte  nanoassemblies  having  tailored 
properties  and  for  development  of  greatly  improved  polyelectrolyte  separation  processes. 
Beyond  its  intrinsic  scientific  importance,  this  project  will  significantly  aid  the  Army  in 
advancing  its  missions  of  developing  new  materials  and  materials  processes,  and 
establishing  fundamental  relationships  between  processing  variables  (here,  fluid  flows, 
electric  fields,  and  surface  patterning)  and  material  behavior.  This  project  also  builds  a 
solid  theoretical  underpinning  for  the  technologies  mentioned  above  (e.g.,  nano  fluidics 
and  microfluidics),  which  are  of  great  interest  to  the  Army  and  promise  reduced  costs, 
smaller  weight  and  size,  and  higher  reliability  over  existing  technologies. 

Brownian  Dynamics  Simulations 

The  principal  tool  used  in  this  project  is  the  Brownian  dynamics  simulation  method. 
Brownian  dynamics  (BD)  is  a  coarse-grained  simulation  method  that  is  capable  of 
describing  polymer  dynamics  in  non-equilibrium  situations  over  relatively  long  time 
scales  (~  1  s).  In  BD,  a  polymer  molecule  is  modeled  as  a  series  of  beads  connected  by 
rigid  rods  (or  flexible  springs),  while  the  solvent  in  which  the  polymer  is  immersed  is 
taken  to  be  a  thermal  bath  that  exerts  random  forces  on  the  beads.  The  polymer  dynamics 
are  followed  by  integrating  equations  of  motion  for  each  bead,  where  the  active  forces  are 
those  due  to  viscous  drag,  the  connectors  between  the  beads,  Brownian  motion, 
externally  applied  flow  and  electric  fields,  and  inter-  and  intra-molecular  interactions. 

The  expressions  for  each  of  these  forces  can  be  related  to  parameters  that  can  be 
measured  experimentally  or  calculated  from  finer-grained  simulation  methods  such  as 
molecular  dynamics.  BD  simulations  are  carried  out  by  performing  a  set  of  runs  and  then 
averaging  over  the  ensemble  to  obtain  properties  of  interest  such  as  the  conformations  of 
molecules  adsorbed  to  a  surface  and  the  velocity  of  a  charged  polymer  in  an  electric  field. 
Although  BD  uses  a  coarse-grained  description  of  polymer  molecules,  it  allows  for 
simulations  over  time  scales  much  longer  than  those  that  can  be  accessed  by  molecular 
dynamics  simulations.  Because  many  phenomena  of  interest  in  polyelectrolyte 
adsorption,  self-assembly,  and  transport  occur  on  such  long  time  scales,  BD  is  an  ideal 
tool  for  studying  them. 


Below,  we  summarize  our  progress  and  accomplishments  on  several  research  problems 
related  to  the  above  project. 


Summary  of  Most  Important  Results 

Polvelectrolyte  Adsorption  in  Shear  Flow 

In  order  to  understand  polyelectrolyte  adsorption  at  a  fundamental  level,  the  logical  place 
to  start  is  the  adsorption  of  single  molecules.  We  performed  the  first  molecular 
simulations  of  polyelectrolyte  adsorption  in  shear  flow  [2],  and  were  the  first  to 
incorporate  hydrodynamic  interaction  (HI)  into  such  simulations  [9,1 1].  The  key  insight 
that  comes  out  of  these  simulations  is  the  discovery  of  a  competition  between  polymer- 
surface  attraction,  which  draws  the  polymer  to  the  adsorbing  surface,  and  polymer- 
surface  HI,  which  drives  the  polymer  away  from  the  adsorbing  surface.  By  combining 
the  simulation  results  along  with  a  kinetic  theory  [7],  we  were  able  to  derive  a  scaling  law 
relating  the  critical  shear  rate  needed  for  desorption  to  various  problem  parameters. 

Our  results  demonstrate  how  shear  flow  influences  the  adsorption  of  polyelectrolyte 
molecules,  and  show  the  critical  role  that  HI  plays.  The  simulation  method  and  results  of 
this  work  are  expected  to  be  useful  for  applications  in  which  polyelectrolyte  adsorption 
occurs  in  flowing  fluids,  such  as  biochemical  sensors,  nanofluidics  and  microfluidics,  and 
biocompatible  coatings  and  materials.  Our  efforts  will  lay  the  basis  for  future  work  in 
which  we  simulate  the  adsorption  and  self-assembly  of  multiple  polyelectrolyte 
molecules  onto  a  patterned  surface  in  the  presence  of  fluid  flow.  As  most  experiments 
involve  multiple  molecules,  we  expect  that  we  will  be  able  to  compare  predictions  from 
our  simulations  with  extensive  experimental  data  already  in  the  literature.  The 
fundamental  studies  on  single -molecule  adsorption  we  have  conducted  will  ultimately 
yield  insight  into  how  to  design  a  flow  field  to  achieve  a  polyelectrolyte  nanoassembly 
with  tailored  properties. 

Polvelectrolyte  Adsorption  onto  Patterned  Surfaces 

Patterning  of  surface  charge  and  topography  is  a  powerful  way  to  influence  the 
conformations  of  adsorbed  polymers.  We  performed  the  most  comprehensive  molecular 
simulations  to  date  on  polyelectrolyte  adsorption  onto  surfaces  decorated  with  periodic 
arrays  of  charged  patches  [5].  We  also  perfonned  the  first  molecular  simulations  of 
polyelectrolyte  adsorption  onto  topographically  patterned  surfaces  [8].  The  key  insight 
that  comes  out  of  these  studies  is  the  discovery  that  the  spatial  structure  of  the  electric 
field  near  the  surface  can  have  a  significant  impact  on  polymer  conformations.  This  is  in 
addition  to  any  confinement  effect  due  to  the  geometry  of  the  charge  or  topographic 
patterns. 

Our  results  demonstrate  how  polyelectrolyte  conformations  can  be  controlled  by  the 
design  of  surface  charge  and  topography .  Current  work  is  focused  on  how  dynamic 
properties,  such  as  diffusion  and  rheology,  are  influenced  by  surface  charge  and 
topography  [13].  We  expect  that  our  studies  will  be  useful  for  applications  where  single 
chain  adsorption  is  important  such  as  sensors  and  biological  processes,  and  also  those 
where  multiple  chains  may  adsorb,  such  as  polyelectrolyte  multilayer  formation.  The 


fundamental  studies  on  single -molecule  adsorption  we  have  conducted  will  ultimately 
yield  insight  into  how  to  design  a  surface  to  achieve  a  polyelectrolyte  nanoassembly  with 
tailored  properties. 

Dendrimer  and  Hyperbranched  Polymer  Adsorption 

Whereas  the  simulations  described  above  considered  only  linear  polymers,  polymers  can 
have  more  complex  architectures  and  these  can  be  useful  for  applications  such  as  surface 
functionalization,  biochemical  sensing,  and  drug  delivery,  all  of  which  are  of  relevance  to 
the  Army.  We  performed  the  first  BD  simulations  of  dendrimer  adsorption,  considering 
both  electrostatic  [6]  and  non-electrostatic  interactions  [12].  Prior  simulations  of  this 
phenomenon  had  been  performed  with  either  finer-grained  simulation  methods 
(molecular  dynamics)  or  equilibrium  methods  (Monte  Carlo).  The  coarse-grained  nature 
of  our  simulations  makes  them  much  less  computationally  intensive  than  molecular 
dynamics  simulations,  and  this  enabled  us  to  perform  the  first  comprehensive 
investigation  of  how  the  distribution  of  charge  and  hydrophobicity  within  dendrimers 
influences  adsorption  behavior.  The  dynamic  nature  of  our  simulations  allows  us  to 
probe  non-equilibrium  behavior  in  a  rigorous  way,  overcoming  a  key  limitation  of  Monte 
Carlo  simulations. 

Our  results  are  expected  to  be  helpful  in  providing  qualitative  guidance  for  dendrimer 
design  in  various  applications,  especially  since  it  is  extremely  cumbersome  to  address  the 
large  parameter  space  involved  through  experiments  alone.  We  are  currently 
incorporating  the  effects  of  fluid  flow,  and  then  plan  to  consider  more  complex 
hyperbranched  polymers  (HBPs)  and  their  self-assembly  with  nanoparticles  for  surface 
functionalization  applications.  These  single-molecule  studies  will  serve  as  a  prelude  to 
simulations  in  which  we  examine  how  multiple  HBPs  self-assemble  with  nanoparticles, 
and  assembly  behavior  near  surfaces  and  in  fluid  flows.  Ultimately,  we  expect  that  the 
results  of  our  work  will  help  guide  experimentalists  to  combinations  of  nanoparticles  and 
HBPs  which  are  promising  for  surface  functionalization. 

Polvelectrolvte  Electrophoresis 

Polyelectrolyte  separations  are  a  crucial  component  in  biosensing,  and  are  significantly 
influenced  by  the  nature  of  the  interactions  between  polyelectrolyte  molecules  and  the 
surrounding  surfaces,  which  may  be  patterned  with  charge  and  topography.  We 
performed  two  BD  studies  of  polyelectrolyte  electrophoresis  through  narrow 
constrictions  [1,3].  Whereas  previous  studies  examined  transport  in  limiting  cases  where 
only  one  time  scale  is  dominant,  the  major  advance  of  our  work  is  the  elucidation  of  how 
multiple  time  scales  can  interact  with  each  other.  The  results  help  explain  existing 
experimental  studies  on  DNA  transport  through  entropic  traps  and  nanopores.  Future 
work  will  focus  on  incorporating  HI  into  the  simulations,  examining  how  chemical 
patterning  of  the  surfaces  can  be  used  to  further  control  transport,  and  the  effect  of 
multiple  polyelectrolyte  molecules.  It  is  important  to  emphasize  that  many  of  the  tools 
and  ideas  developed  in  studying  polyelectrolyte  adsorption  and  self-assembly  can  be 
directly  applied  to  study  electrophoresis  problems.  We  expect  that  the  fundamental 
understanding  gained  from  our  studies  will  provide  insight  into  how  to  design  surfaces  to 


optimize  polyelectrolyte  separations,  which  in  turn  will  advance  the  development  of 
biosensors. 

Continuum  Modeling  of  Adsorption-related  Phenomena 

In  addition  to  the  molecular-level  simulations  described  above,  we  also  developed 
continuum-level  models  of  various  adsorption-related  phenomena.  Such  models  are 
valuable  for  several  reasons.  First,  they  can  sometimes  be  fonnulated  and  solved  quickly 
to  obtain  insight  which  can  then  be  used  to  perform  computationally  intense  molecular- 
level  simulations  more  efficiently.  Second,  if  the  insights  gained  from  molecular-level 
simulations  can  be  distilled  into  a  continuum-level  model,  then  the  model  can  serve  as  a 
convenient  tool  for  prediction,  design,  and  optimization  in  engineering  applications. 
Third,  by  studying  problems  from  both  molecular-  and  continuum-level  perspectives, 
knowledge  can  be  gained  about  how  to  develop  multiscale  simulation  tools  that  bridge 
the  gap  between  small  and  large  scales. 

During  this  work,  we  developed  several  continuum-level  models: 

•  The  first  kinetic  theory  of  polyelectrolyte  adsorption  in  shear  flow  [7].  The  insight 
gained  was  instrumental  in  developing,  validating,  and  interpreting  the  BD  simulations 
discussed  above,  which  can  address  situations  too  complex  for  the  kinetic  theory  to 
handle. 

•  The  first  mean-field  theory  to  study  polyelectrolyte  adsorption  onto  surfaces  patterned 
with  charge  and  topography  [10].  The  predictions  of  the  theory  provided  considerable 
additional  physical  insight  into  the  BD  simulation  results  discussed  above,  which  in  turn 
can  be  used  to  reveal  a  number  of  phenomena  beyond  the  scope  of  the  theory. 

•  The  first  free-energy  model  of  polyelectrolyte  adsorption  onto  topographically 
patterned  surfaces  [8].  The  model  enabled  us  to  rationalize  the  BD  simulation  results 
discussed  above,  and  may  be  useful  to  other  researchers  in  the  field. 

•  The  first  linear  stability  analysis  of  a  power-law  fluid  flowing  past  a  soft  polymer  gel 
[4].  The  analysis  was  motivated  by  the  experimental  observation  that  polyelectrolyte 
multilayers  have  rheological  properties  very  similar  to  soft  polymer  gels.  Since  fluids 
flowing  past  soft  polymer  gels  can  become  unstable  and  lead  to  degradation  of  the  gel, 
multilayers  may  be  susceptible  to  a  similar  degradation.  Our  analysis  showed  that 
instability  is  unlikely  to  occur  for  the  range  of  problem  parameters  relevant  to  multilayer 
applications.  Thus,  if  multilayer  degradation  in  a  fluid  flow  is  observed  experimentally, 
it  will  likely  be  because  of  a  completely  different  physical  mechanism. 

Overall,  these  models  complemented  the  knowledge  obtained  from  our  molecular 
simulations,  deepened  our  physical  understanding  of  the  simulation  results,  and 
demonstrated  how  molecular-  and  continuum-level  approaches  can  be  combined  in  a 
synergistic  way. 
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